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Abstract 

The purpose of this work is to demonstrate how an arbitrarily chosen background of the Universe 
can be made a solution of a simple geometric sigma model. Geometric sigma models are purely 
geometric theories in which spacetime coordinates are seen as scalar fields coupled to gravity. 
Although they look like ordinary sigma models, they have the peculiarity that their complete 
matter content can be gauged away. The remaining geometric theory possesses a background 
solution that is predefined in the process of constructing the theory. The fact that background 
configuration is specified in advance is another peculiarity of geometric sigma models. In this paper, 
I construct geometric sigma models based on different background geometries of the Universe. 
Whatever background geometry is chosen, the dynamics of its small perturbations is shown to 
posses a generic classical stability. This way, any freely chosen background metric is made a stable 
solution of a simple model. Three particular models of the Universe are considered as examples of 
how this is done in practice. 
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I. INTRODUCTION 


The latest astronomical observations have given a substantial boost to the development 


of modern cosmology 


B 


10l |. In particular, the accelerating expansion of the Universe has 


drawn much attention. The early time acceleration is widely known as inflation, while the 


late time acceleration is usually referred to as the epoch of dark energy lll-ll4j|. Presently, 


the ACDM model, in which the cosmological constant A plays the role of dark energy, is 
accepted as a standard cosmological model. There is an extensive literature on other forms 


of dark energy, too 


15M25l|. All in all, the number of dark energy models that can be 


found in literature is enormous. The same holds for the inflationary models that have been 
constructed over the years. 

In this paper, I shall describe the procedure which associates an action functional with an 
arbitrarily chosen background geometry of the Universe. Precisely, any desirable geometry 
of the Universe is made a solution of a particular geometric sigma model. Geometric sigma 
models are theories that possess two distinctive properties. First, their complete matter 
content can be gauged away. Second, any predefined geomet^_can be made a solution of a 
particular model. These models have first been proposed in 2^ in the context of fermionic 
excitations of flat geometry. Here, I use them for modeling the dynamics of the Universe. 
To be more accurate, only geometry and dark energy are considered in this approach. The 
inclusion of ordinary matter is discussed separately. 

The results obtained in this paper are summarized as follows. First, a class of purely 
geometric dark energy models has been constructed. Every particular model is defined 
as a geometric sigma model associated with a spatially flat, homogeneous and isotropic 
geometry. This way, an arbitrarily chosen geometry of this kind becomes a background 
solution of a particular geometric sigma model. Ultimately, one is provided with the class 
of dark energy models parametrized by their background geometries. The inflation and the 
late time acceleration have purely geometric origin. It is important to emphasize that, while 
the background metric can be chosen arbitrarily, the physics of its small perturbations can 
not. In fact, the background metrics just parametrize geometric sigma models, very much 
the same as inflaton potentials parametrize the inflationary models. 

The second result concerns the linear stability of the background solution in geometric 
sigma models. It has been proven true for almost all background geometries. Precisely, the 
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stability is guaranteed up to the existence of critical moments, where the perturbations may 
diverge. There, however, the linear analysis fails, and should be corrected by the inclusion 
of interaction terms. 

Finally, I have analyzed geometric sigma models coupled to ordinary matter. It has been 
shown that matter helds do not compromise the vacuum stability established earlier. In the 
case of minimal coupling to the metric, the linear stability of matter itself has been proven. 

The layout of the paper is as follows. In Sec. [TTl the construction of geometric sigma 
models, as dehned in [2^, is recapitulated and subsequently applied to spatially flat, homo¬ 
geneous and isotropic geometries. As a result, a class of action functionals of the Universe 
is obtained. Each of these action functionals possesses a nontrivial background solution 
that describes the background geometry of a particular Universe. In Sec. IIIIl the dynamics 
of small perturbations of these nontrivial backgrounds is examined. In Sec. IIVI the back¬ 
ground solutions are proven stable for almost all spatially flat, homogeneous and isotropic 
geometries. In Sec. El geometric sigma modes are coupled to ordinary matter. It is shown 
that matter helds preserve the results obtained in the absence of matter. In Sec. I VII the 
examples of inhationary and bouncing Universes are used to demonstrate how geometric 
sigma models are constructed in practice. Sec. IVIII is devoted to concluding remarks. 

My conventions are as follows. The indexes /x, u, ... and i, j, ... from the middle of 
the alphabet take values 0,1,2,3. The indexes a, f3, ... and a, b, ... from the beginning of 
the alphabet take values 1, 2, 3. The spacetime coordinates are denoted by x^, the ordinary 
diherentiation uses comma = 9^X), and the covariant diherentiation uses semicolon 

(X^ = V^X). The repeated indexes denote summation: Xaa = Xu -|- X 22 + X 33 . The 
signature of the 4-metric is (—, -|-, -|-, -t-), and the curvature tensor is dehned as R^uXp = 
d\Tf^^p — dpT^^^x + 


II. GEOMETRIC SIGMA MODELS 


The construction of a geometric sigma model begins with specifying a spacetime metric. 
I shall denote it with The metric g^^J is freely chosen, and the coordinates are 

fully hxed. As a consequence, the functions gj^u {x) are completely determined. In the next 
step, the corresponding Ricci tensor R^p^ix) is calculated, and the following Einstein like 
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equation is postulated: 


= ( 1 ) 

Obviously, the metric gjfj is a solution of the equation flTJ. Its non-zero right hand side 

defines matter content of the theory. The equation ([T]) is an example of how an arbitrarily 
chosen metric can be made a solution of a simple model. 

The equation ([1]) obviously lacks general covariance. To covariantize it, I introduce a new 
set of coordinates (jf = 0*(x). In terms of these new coordinates, the equation ([1]) takes the 
form 

( 2 ) 

where the functions are dehned through 

= (3) 


In other words, the ten functions are obtained by replacing x with 0 in ten components 

of the Ricci tensor R^^}{x). The equation ([2]) is generally covariant once the new coordinates 
0* are seen as scalar functions of the old coordinates x^. If the new coordinates are chosen 
to coincide with the old ones, 0*(a:) = the covariant equation ([ 2 ]) is brought back to 

the non-covariant form ([1]). In what follows, the shorthand notation = x* will be used. 

The equation ([2]) has the form of the Einstein’s equation in which four scalar helds <fi{x) 
of some nonlinear sigma model are coupled to gravity. The ’’matter held equations” are 
obtained by utilizing the Bianchi identities — ^g^'^R).i, = 0. If the condition det 0*^ 7 ^ 0 
is fulhlled, one obtains 


-I- i 



an it dHtA 

00 * 00 -? / 



0 . 


(4) 


The equation (jTj) is not an independent equation, as it follows from ([2]) and the Bianchi 
identities. It is straightforward to verify that the equations ([2]) and (jl]) can be derived from 
the action functional 

4 = ^ / [R - i/,, (0)0>^’^] . (5) 

Here, the target metric is not an independent coefficient of the model. Instead, it is 

constructed out of the background metric through its dehning relation ([2]). This way, 
an action functional is associated with every freely chosen background metric. This action 
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functional describes a nonlinear sigma model coupled to gravity, and possesses the nontrivial 
(let me call it vacuum) solution 

(f>" = x\ g^u = ■ ( 6 ) 

Indeed, the equation ([2]) with the target metric ([3]) is trivially satished if the scalars 0* and 
the metric g^i, are given by ([ 6 ]). Being a direct consequence of ([2]), so is the equation (jl]). 

The physics of small perturbations of the vacuum ([ 6 ]) does not violate the condition 
det 0*^ 7 ^ 0, which enables one to interpret the scalars 0* as spacetime coordinates. If this is 
the case, one is allowed to fix the gauge 0 *(a:) = a:*, which brings us back to the geometric 
equation ([ 1 ]). One should have in mind, however, that the action ([5]) has a non-geometric 
sector, too. It is characterized by det 0(^ = 0, and includes trivial vacuum solutions such as 
0* = const. In what follows, I shall restrict my considerations to purely geometric dynamics 
of small perturbations of the vacuum ([ 6 ]). 

Before I continue, let me note that the equation ([T]) is not the unique geometric equation 
that allows the solution = gj^J. A simple generalization of this equation can be obtained 
by adding terms proportional to gf^i, — gj^!. The simplest choice is the equation 

V = + Iv'(x) . (7) 

It defines a class of geometric theories parametrized by metrics gl)^J, and potentials V. The 
covariantization of the non-covariant equation ([7]) ultimately leads to the action functional 

4 = ^ / ‘1‘xs/^ [fi - - V(^)] , (8) 

where the target metric Fij{(f)) is dehned by 

Fij{x) = {x)gff (x). (9) 

The class of theories dehned by ([8]) possesses the vacuum solution ([6]) for any choice of 
the potential ld(0). The physics of small perturbations of this vacuum allows the gauge 
condition 0* = x\ which brigs us back to the geometric equation ([7]). 

In what follows, I shall associate a geometric sigma model with a given vacuum geometry 
of the Universe. Let me choose a spatially hat, homogeneous and isotropic metric gl^J, 
dehned by 

-I- a‘^{t) [dx^ + dy^ + dz^) . (10) 
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The corresponding Ricci tensor is calculated straightforwardly. One obtains 

= -3 - , Rtl = 0 , = {ad + 2d^) , 


^-00 

where ’’dot” denotes time derivative. Now, I am ready to construct the target metric Fij{(j)). 
In this paper, I choose the simplest model allowed by ([9]). It is obtained by noticing that 
the target metric Fij can signihcantly be simplihed by a proper choice of the potential V. 
Indeed, if the potential is chosen to have the form 


V(t)^2l2^ + - 


( 11 ) 

the component Tqo remains the only nonzero component of the target metric. Precisely, one 
obtains 

, , 

F{t = 2 - - - , 12 

\a^ a J 

where the identihcation Fqo = F is introduced for convenience. As a consequence, 0° is the 
only scalar held that enters the action functional ([8]). Let me simplify the notation by using 
the identihcation (jp = 0. The action ([8]) then reduces to 

1 






(13) 


It governs the dynamics of gravity coupled to a scalar held, and possesses the vacuum 
solution 

0 = ^, g^2v = ■ (14) 

The precise form of the target metric R(0) and the potential R(0) is determined once the 
function a{t) is specihed. It should be noted that it is only the background geometry gjpj 
that is freely chosen. The dynamics of metric perturbations is governed by the corresponding 
action functional. The class of action functionals flT3l) represents a collection of dark energy 
models parametrized by the scale factors a(t). 


Similar attempts to derive a dark energy moc 
in literature. Take, for example, references 27 


el out of the given scale factor already exist 
. There, scalar held dark energy models 


have been constructed to mimic holographic dark energy. Every particular model is built 
with a separate ehort to solve a particular problem. The procedure described in this paper, 
however, is the hrst systematic approach of the kind. It gives a precise prescription of how to 
construct the target metric and the potential of a stable dark energy model. As we shall see 
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later, the generic classical stability is guaranteed for nearly any chosen background. Some 
specihc models are considered in Sec. IVTl 

It should be noted that no ordinary matter has been considered so far. Luckily, the 
inclusion of matter helds does not compromise the basic predictions of dark energy models 
([12]). This will be demonstrated later in Sec. jVj Besides, the prevailing form of matter in 
the Universe is believed to be the dark energy. Thus, the class of dark energy models (fT3|) 
can roughly be interpreted as zero approximation of more realistic cosmologies. 

The standard physical requirements that ensure the absence of ghosts and tachions re¬ 
strain the target metric F{(j)) to be positively dehnite, and the potential U(0) to be bounded 
from below. These restrictions, however, refer to trivial vacuums. Precisely, the positive 
dehniteness of and the fact that (j) = (po is a. minimum of the potential V{(j)) ensure sta¬ 
bility of the vacuum (p = (po, In this paper, however, the vacuum of interest is the 

nontrivial vacuum fUTj) . Its stability is not guaranteed by the above physical requirements, 
and I am led to check it by direct calculation. 

III. DYNAMICS OF SMALL PERTURBATIONS 

In this section, I shall examine the dynamics of small perturbations of the vacuum ffTT|) . 
as governed by the action functional (1T3|) . The inhnitesimal change of coordinates 

leaves this action invariant, and allows the gauge hxing cp = t. In this gauge, 
the matter held equation is identically satished, and we are left with the gravitational held 
equation ([7]). The residual diheomorphisms are dehned by the constraint = 0. 

The only variable in the gauge hxed theory is the metric perturbation h^y, dehned by 

9 iiu — 9I1U F h^y. 

With respect to the residual diheomorphisms, it transforms as 

<^0^00 = 0 , 

® (^o,/3 T ^13,a) y 

where 5o is the form variation, and It is seen that hoa can also be gauged away. 

The gauge condition 

hoa 0 
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restrains the gauge parameters to be functions of spatial coordinates, only. Precisely, the 
residual gauge parameters are dehned by 

e° = o, r = o. 


In what follows, I shall demonstrate how the residual gauge symmetry can further be hxed. 

Let me hrst linearize the held equations ([7]). After cumbersome, but straightforward, 
calculation one obtains 


do I haa — 2— j + 3(2(2 /iQO A 2 (ah + 2(2^) hoo + hoo^aa — 0 : 


do 






2— hoo,a — 0 , 
a 


■■ a ■ a 1 

2 hcn^ “ 1 “ 2 (^ 0 : 7 ,^/37,7a h, 

_ CL CL CL 

CL ' CL^ 

— h,y,y — 2 —— /i™ OjOj hoo 2 (cLQj 2(2^) hoo 

_a a^ ^ 

It is immediately seen that the equation fllSbj) implies 

1 


'al3,'y'y 

hal3 + hoo,0/3 = 0 


(15a) 

(15b) 

(15c) 


do 


(ho7,7/3 hjS'y, 


70 J 


= 0 , 


which tells us that the expression in square brackets does not depend on time. As a conse¬ 
quence, this expression can be gauged away. Indeed, its transformation law reads 

1 


(^07,7/3 


^a,'y'y/3 'C/3,' 


,770 • 


Both, the expression in square brackets and the residual gauge parameters Co are functions 
of spatial coordinates alone. This allows the gauge hxing 


^■ 07 , 7/3 0. (16) 

In what follows, I shall simplify the analysis by the assumption that metric perturbations are 
spatially localized. This means that the perturbation h^^ is assumed to decrease sufficiently 
fast in spatial inhnity. With this assumption, many expressions are simplihed. For exam¬ 
ple, the equation X^a = 0 has the general solution X = X{t), but the adopted boundary 
conditions imply X = 0. Similarly, the equation Xq,„ = 0 has the unique solution X = 0. 
With this in mind, one easily determines the residual symmetry after the gauge condition 
([TB]) has been imposed. It is dehned by 


C« = e,' 
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where the new parameter e is an arbitrary function of spatial coordinates. 

Let me now extract the divergence free parts of the variable hap- To this end, I use the 
decomposition 

hap — hap T ha^p T hp a T h ap i 

where hap and ha are, by dehnition, divergence free {hap,p = ha^a = 0). In what follows, 
all the expressions will be rewritten in terms of the new variables hap, ha and h. Let me 
start with the gauge condition flTHl) . With the help of the adopted boundary conditions, it 
is straightforward to verify that flTB]) becomes 

ha = 0. 


Now, I am ready to rewrite the equations flT^ in terms of the remaining variables hgo, 
hap and h. Let me hrst consider the scalar sector. As it turns out, there are only three 
independent scalar equations. The hrst two are the constraint equations 


ad 

o ^OLOL 

a 

+ - 

—h 

= 2a^ 

4^ 


_a^ 

,0 ® 

_a^ 

,/3/3 

_a^ 


,0/3/3 


1 


hr 


J ,0 


+ 2— hoo — 0 . 
a 


(17a) 


(17b) 


The equation fllTbl) follows from fllSbl) . and fll7al) is a linear combination of the trace and 
divergence of (libel) . It is seen that hoo is fully determined by other variables. Thus, it 
carries no degrees of freedom. On the other hand, the variable h/a^ is determined up to a 
free function of the spatial coordinates. This freedom, however, can readily be gauged away. 
Indeed, the variable h/a? transforms as 

hn ( —h 1 = —2e 


with respect to the residual symmetry of the model. As haa/o? is gauge invariant, the 
residual parameter e = e{x) is exactly what one needs to hx the free integration function 
of fll7ap . In summary, the constraint equations (fT7|) tell us that neither hoo nor h/a? carry 
physical degrees of freedom. 

The third scalar equation is the equation fllbap . or equivalently, the divergence of (libel) . 
When its coefficients are expressed in terms of the Hubble parameter H = a/a, it takes the 
form 



’ 1 r 

/ • 0 -0 -A 

■ 1 ~ 

1 

■ 1 ~ 

HH 

_ h 

0 

+ -2H^ + HHj 

,00 ^ ^ 

0 

- ^hh 

,0 « 

0 ^OlOL 


= 0 . 


(18) 


,/3/3 
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The equation ffTSl) governs the dynamics of the unique scalar mode of the model. 

Now, I am left with the traceless, divergence free part of the equation fllScp . which governs 
the dynamics of the traceless, divergence free part of ha^. The latter is dehned by 


1 ~ 1 

ha 3 — ^a /3 ^77*^0/? T ^ 


77 / 5 


where stands for the inverse of the Laplacian A = (The existence of A“^ 

is guaranteed by the adopted boundary conditions, which state that the perturbations 
decrease sufficiently fast in spatial inhnity.) Then, the tensor part of the equation fllScp 
takes the simple form 


■ 1 ^ 

a 

■ 1 ^ 

1 

■ 1 ^ 


+ 3- 

'^hap 

~Z2 



,00 ® 


,0 « 



= 0 . 


(19) 


J ,77 


Being subject to the constraints haa = = 0, the variable hap carries two physical 

degrees of freedom. 

Before I proceed, let me draw your attention to the fact that the scalar equation flTSll is 
identically satished if if = const. This corresponds to the choice a oc If this is the case, 
the model flT^ reduces to GR with the cosmological term—the model that carries only two 
physical degrees of freedom. In what follows, I shall restrict my considerations to the case 
H 7 ^ const. Then, the scalar equation (ITSll is rewritten as 


2 ^aa 


,00 




:hn 


,0 


2 ^aa 


= 0 


( 20 ) 


-I ,/3/3 


It is seen that the scalar mode of the geometric sigma model flT^ is massless. 

IV. STABILITY ANALYSIS 

In this section, I shall examine the stability of the vacuum solution = 0. Let me start 
with the tensor equation flT^ . First, I introduce the collective variable 


Q — 2 ^ 0/3 • 

Qj 

In terms of Q, the equation fllOp takes the compact form 

Q + 3 — Q l^Q^aa = 0 . 

a 


( 21 ) 
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The function Q{x) is searched for in the form 

Q = Re J d^kq{k,t)e^^'^, (22) 

whereupon the equation fl2TD becomes 

q + 3-q + ^q = 0. (23) 

a 

The stability of the vacuum solution g = 0 is examined by the canonical analysis of the 
equation fl2^ . In the hrst step, I notice that the equation fl23|) is obtained from the La- 
grangian 

C = a {a^q^ — k‘^q‘^) ■ 

Indeed, it is easily verihed that its variation leads to the equation (12^ . The corresponding 
Hamiltonian is straightforwardly calculated to be 

= a + • (24) 

It is seen that the Hamiltonian is positive for all the allowed values of the parameter a{t), and 
all the values of the wave vector k. Its minimum is located at g = p = 0. For Hamiltonians 
with no explicit time dependence, this would imply the stability of the vacuum g = p = 0. 
Indeed, owing to "H = = 0, the physical phase space trajectories coincide with the 

orbits "H = const. These orbits, on the other hand, are closed curves around the vacuum 
g = p = 0. As a consequence, the phase space trajectory which is initially close to the 
vacuum continues to be in the vicinity of the vacuum at all times. 

Unfortunately, the Hamiltonian fl2T)) depends on time through the free parameter a{t). 
As a consequence, 

and the Hamiltonian is not conserved. Still, the stability of the vacuum g = p = 0 is not 
compromised. To see this, note that the orbits T-i = const, remain to be closed curves around 
g = p = 0, only this time they evolve in time as shown in Fig. [1] It is seen from fl2^ that 
the change of the Hamiltonian along the phase space trajectory (g(t),p(t)) is the same as for 
the still point (g,p) = const. Thus, the general solution of fl2^ is the phase space trajectory 
which, at any time, touches the respective orbit "H = const. As these orbits are closed curves 
around g = p = 0, the phase space trajectories initially close to g = p = 0 remain to be close 
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FIG. 1: Time evolution of orbits H = const. 


to g = p = 0 at all times. In other words, the vacuum solution q = p = 0 is stable against 
small perturbations. This holds true for any a{t) and k. 

The stability of the scalar equation fl2U]) is examined analogously. Using the notation 


Q — 9 , 


the equation (j2nil is rewritten as 


Q + ( 3H — 2— + Q — -^Q^aa — 0 , 


and its Fourier transform, as defined by ([22]), becomes 


(26) 


+ = (27) 

(This notation should not be confused with the same notation used in the analysis of tensor 
modes.) It is easily checked that the equation fl27ll follows from the Lagrangian 


H 


jC oc a— {a?q^ - k‘^q^) , 


or equivalently, from the Hamiltonian 


•ay ^ 
n oc a— 


/ 1 



(28) 


The proportionality sign reminds us that C and "H are determined only up to a multiplicative 
constant. If necessary, this freedom can be used to correct the overall sign of the Hamilto¬ 
nian. As opposed to the Hamiltonian fl2T|l . the Hamiltonian (j28jl is not necessarily positive. 
However, it is only the overall sign of 'H that can be negative. This implies that the vacuum 
q = p = 0 can only be a maximum or a minimum, and never a saddle point. As a conse¬ 
quence, the lines of constant 'H are closed curves around g = p = 0. From this point on, 
the stability analysis reduces to that of the tensor mode, leading to the conclusion that the 
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vacuum q = p = 0 is stable against scalar perturbations, too. This holds true in a generic 
time interval, and for a generic choice of a{t). 

Let me note, however, that there can exist critical moments in the evolution of scalar 
perturbations, where the stability may be lost. As seen from fl25]) . these are dehned by iL = 0 
or H = 0. In fact, the stability discussed in this section is proven only up to the presence 
of critical moments. Such critical moments can be found in many cosmological models, as 
will be demonstrated in Sec. lYD Some of them are benign, but some may cause divergent 
behavior. It is important to realize that the construction of geometric sigma models does 
not guarantee the absence of critical moments. As a consequence, not every choice of the 
scale factor a{t) leads to an everywhere regular model. One should have in mind, however, 
that the linear analysis considered in this paper is inapplicable in the vicinity of critical 
moments. Indeed, the perturbations grow big there, and the interaction terms can not be 
neglected. 

Before I close this section, let me mention once more that matter helds have been excluded 
from the above analysis. As a consequence, the proven stability may be compromised. 
Indeed, the target metric F{(j)) has been allowed to take negative values, and the potential 
V{(j)) to be unbounded from below. Then, the conventional coupling to matter helds may 
lead to the appearance of ghosts and tachyons. In this paper, however, I am examining the 
linear stability alone. It will be demonstrated in the next section, that matter helds can not 
appear in the linearized Einstein’s and scalar held equations. Thus, the established stability 
is not threatened. The stability of matter perturbations themselves, on the other hand, is 
studied separately. 

V. MATTER FIELDS 

In this section, I shall examine how the presence of matter helds inhuences the behavior 
of geometric sigma models. My starting point is the action 

I = Ig + Im 1 (29) 

where Ig is the geometric action flT^ . and Im stands for the action of matter helds. Cus¬ 
tomarily, the matter Lagrangian is thought of as the Lagrangian of the standard model of 
elementary particles, minimally coupled to gravity, and possibly, to the inhaton held 0. In 
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what follows, the matter fields will collectively be denoted by hi. 

Let me hrst consider the simplest case characterized by the absence of direct matter- 
inflaton couplings. Instead, the matter helds are minimally coupled to the metric alone. 
This choice is justihed by the fact that minimal coupling to the metric already contains a 
simple coupling to the inflaton itself. Indeed, it will shortly be shown that matter helds 
do not compromise the gauge hxing procedure of the preceding sections. This procedure 
leaves us with three dynamical metric components, one of which is a reminiscent of the 
original inhaton held. Of course, one can always change this simple scenario by employing 
direct couplings to the inhaton held. It will be shown later that most of the direct inhaton 
couplings preserve the results of the preceding sections. 


A. Vacuum solution 


In this subsection, I shall examine how the sigma model vacuum flTTll is inhuenced by the 
presence of matter helds. Let me start with the analysis of the matter held equations 




0 . 


(30) 


Owing to the minimal coupling to the metric, these equations possess vacuum solution that 
coincides with the well known vacuum of the standard model of elementary particles. Indeed, 
the standard model equations are trivially satished when matter helds take proper constant 
values 

hi - , 


and so are the equations (jHilil . This holds true for any value of the metric that appears in 
the held equations. In particular, the vacuum value of the stress-energy tensor 


rjijJLU _ 

m — 


2 5Im 

\/^ 


is zero, irrespectively of the presence of Formally, 


(31) 


n = ^ TZ = o 


for any With this, the inhaton and Einsten’s equations take the form of the sigma 
model equations considered in the preceding sections. Indeed, owing to the absence of the 
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matter-inflaton coupling, the inflaton equation is the same as that of the sigma model, 


51 

6 ^ 


6Ig 

6(f) 


= 0 . 


Einstein’s equations, on the other hand, reduce to the sigma model equations once the 
matter helds take their vacuum values. 




It has already been shown in Sec. [TTl that geometric sigma model possesses the vacuum 
solution (HI. As a consequence, the held equations that follow from the action (|2^ are 
satished by 

= ffo , 0 = t, gi,v = glfff} ■ (32) 

This is the vacuum, or shall we say, the background solution of the model fl29l) . It is seen 
that the presence of matter helds does not compromise the sigma model vacuum of the 
preceding sections. This result is summarized in the sentence 


• matter fields do not violate the sigma model vacuum. 

In what follows, I shall examine the stability of the vacuum fl32l) against its small pertur¬ 
bations. Owing to their smallness, the dynamics of vacuum perturbations is governed by 
the linearized held equations. Thus, it is the linear stability of the model fl2^ that will be 
examined in the next subsection. 


B. Stability analysis 

The linear stability of the vacuum fl32l) is examined by inspecting the linearized held 
equations of the action fl29D . It is immediately seen that, after linearization, the inhaton 
and Einstein’s equations reduce to those of the geometric sigma model of the preceding 
sections. Indeed, the stress-energy tensor Tfffi, being at least quadratic in perturbations 
of matter helds, does not appear on the r.h.s. of the linearized Einstein’s equations. At 
the same time, the inhaton does not couple to matter helds, at all. Hence, the linearized 
inhaton and metric equations of motion remain unchanged by the inclusion of matter. They 
are diheomorphism invariant, so that the complete gauge hxing procedure of Sec. Illll is still 


15 


valid. In this gauge, the tensor and scalar equations (IT^ and (EO]) are exactly what one 
obtains from the linearized equations 6I/6(f) = 0 and 51 /5gfj,u = 0. 

The generic linear stability of the sigma model vacuum has already been established in 
the preceding section. What remains to be shown is the linear stability of the matter held 
equations. To this end, notice that the minimal coupling to matter helds implies the validity 
of the principle of equivalence. Indeed, the matter helds of the held equations SI/SQ = 0 
are coupled to at most hrst derivatives of the metric gn^. In a local inertial frame, the 
metric derivatives vanish, and the metric itself becomes the Minkowski metric rj^u- This 
way, the matter held equations take their special relativistic form—the standard model of 
elementary particles in hat spacetime. The latter is known to be stable against perturbations 
of its trivial vacuum. As a consequence, the linearized matter held equations possess the 
stable vacuum = flo- To summarize, the vacuum fl3^ has a generic linear stability against 
perturbations governed by the action (1291) . The dynamics of its geometric part remains the 
same as found in Sec. IIIIl Therefore, 

• the presence of matter fields does not violate the established linear dynamics of geo¬ 
metric sigma models. 


What remains to be found is the dynamics of the linearized matter held equations. There 
are three types of matter helds that appear in the action: fermion helds, gauge helds and 
the Higgs. In the linearized theory, fermion helds are governed by the Dirac equation, gauge 
helds obey Maxwell equations, and the Higss is subject to the Klein-Gordon equation. All 
of these are minimally coupled to the external curved background. 

Scalar fields. The scalar (p, minimally coupled to the external metric gjf,}, obeys the 
Klein-Gordon equation (□ —= 0. The scalar held mass is denoted by m, and □ stands 
for the covariant d’Alembertian of the vacuum metric gjf,}. With gjf^ of the form (1T0|) . the 
Klein-Gordon equation becomes 

ip + 3-(p (m^ —= 0 ) 

a \ J 

where A = Now, one can use the Fourier decomposition of the preceding 

section to rewrite the above equation in the form 
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This equation differs from the scalar perturbation equation fl27D not only by the presence 
of mass, but also by the different friction coefficient. During the inflationary phase, when 
Hubble parameter is approximately constant, the friction coefficient of fl27p is, in most cases, 
close to that of fl5^ . In some cases, however, the scalar modes of the geometric sigma model 
may have significantly higher friction, as will be demonstrated in the next section. In such 
situations, the rapid expansion of the Universe makes the inflaton decay much faster than 
scalars of the matter Lagrangian. 

Gauge fields. In linear approximation, gauge fields obey the equation V= 0, 
where is the linearized gauge held strength, and stands for the 

covariant derivative. To simplify calculations, I shall work in the Coulomb gauge. Then, the 
time component Aq is constrained to be zero, while the spatial components A^ satisfy 

Aa H— Ao, - -AAa = 0 , do,Aa = 0 . 

a 

In terms of their Fourier components Qa, the held equations take the form 

a 

Qa + = 0 , kaQa = 0 . (34) 

a 

It is seen that the friction coefficient of (l34)l is three times smaller than that of (IMjl . Thus, 
in the expanding Universe, the gauge helds outlast the excitations of the gravitational held. 

Dirac field. The evolution of Dirac held minimally coupled to gravity is given by the 
equation 

k - m)'ip = 0 , (35) 

where 7 ^ are Dirac gamma matrices ({ 7 *, 7 -^} = — 277 ®-^), and stands for the covariant 
derivative 

The gravitational variables that enter this equation are the tetrad and the spin connec¬ 
tion The tetrad is the inverse of while o;®-^^ are the components of in the 

basis of Lorentz generators | [ 77 ,7j], 

1 o' 

= -UJ \(Tij . 

The metric and the connection which are used in this paper, are related to 

and cu®-^^ through the equations 


^ A uj n T ni/^ \ 0. 


(36) 
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The first equation is the very dehnition of the orthonormal tetrad, while the second represents 
the metricity condition. It is seen that, given the metric the tetrad is determined only 
up to the local Lorentz rotations. These, however, are a symmetry of the Dirac equation 
035]), and can be gauge hxed. Thus, starting with the background metric , the simplest 
solution for the background tetrad is found to have the diagonal form 

e°o = l, e\ = e ^2 = e ^3 = a. 


Then, the second equation fl36|l yields the background value of the spin connection. The 
only non-zero components turn out to be 

, _ _ Ob _ -cb 

^ a — ^ a — * 


With the known background values of the tetrad and spin connection, the Dirac equation 
is rewritten as 

-0 + -7°7"V’,a + 7 ( 3 - + 2im7° ] •0 = 0. 
a 2 \ a / 


Its Fourier expansion then yields 


^+2 


3- -|- 2 * 7 ° ( m -\—7 ■ k 


a 


a 


q = 0, 


(37) 


where q{k,t) are dehned by 


'4){x) = / d^k q{k, t) e 


ik-x 


Using the notation 


u \ / «+ \ / W 

q = \ I, u = \ I, V = 


VI \ “- / \ 

the 4-component equation (l37ll is rewritten as a pair of 2-component equations. The hrst is 
the constraint equation 


V = 


m 


u -\ -( 3a — 2id ■ k] u 

2a 


which tells us that v carries no degrees of freedom. The second is the dynamical equation 


il -I- 3Hu + 


— ■ ( —h iHd] — ( 2H + 3H^ 
a \ a / 4 


u = 0. 


If the coordinates are rotated so that k is directed along the 2 ;-axis, the above equation turns 
into a system of two one-component equations 


u± -|- 3Hu± + 


+ ^ i2H + 3H^ 
a \ a A 


= 0. 


(38) 
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It is seen that the friction coefficient coincides with that of the scalar held, whereas the mass 
term is different. In the early stage of inhation, when a is still very small, the inequality 
k ^ aH holds for a wide range of k. In this regime, the mass term is close to that of the 
scalar held. At the end of inhation, however, the scale factor is large, so that most k satisfy 
k -C aH. Then, the mass term takes the form + (3if/2)^. As the inhationary value of H 
is typically much larger than masses of elementary fermions, the equation (l38|) is practically 
independent of m. This means that the production and propagation of elementary fermions 
are the same for all the fermion species. 

C. Matter-inflaton coupling 

So far, the considered matter helds have been assumed to couple to the metric alone. 
It has been shown that the results of the previous sections are not compromised by the 
inclusion of such matter. In this subsection, I shall discuss direct matter-inhaton couplings. 

Let me start with the verihcation of the vacuum solution fl32|) . It is immediately seen 
that interaction terms which are at least quadratic in perturbations of matter helds do not 
compromise the vacuum ([32]). Indeed, their contribution to the held equations disappears 
when 'ip = 'ipQ, thereby making (l3^ a valid vacuum solution. The sigma model vacuum 
oa. on the other hand, satishes the held equations even when the couplings are linear in 
perturbations of matter helds. Thus, 

• typical matter-inflaton couplings preserve the sigma model vacuum. 

The examples of possible matter-inhaton couplings are 

and so on. It is seen that symmetries of the standard model Lagrangian restrict most of 
these interaction terms to be at least quadratic in matter helds. This is fortunate because 
such matter-inhaton couplings preserve the results of the preceding sections. Indeed, the 
linearized Einstein’s and inhaton equations remain the same as those obtained in the absence 
of matter helds. Thus, 

• the analysis of the preceding sections is not compromised by the presence of interaction 
terms which are at least guadratic in matter fields. 
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What does change, however, is the form of the linearized matter held eqnations. As an 
example, let me consider Dirac held with the inhaton conpling of the form /(0)'0'0. Then, 
the linearized eqnations fl38p take the form 


u± + Au± + 

where M{t) = m — f(t), and 


+ + ^H {2C + 3H) 


u± = 0 , 


A = 3H- — , B = H+ — , C = — - — . 

M’ M ' H M 


(39) 


It is seen that the eqnation (l39|) has singular points dehned by M(t) = 0. The easiest way 
to get rid of these is to restrain the couplings to satisfy /(0) < m. Such is, for example, 
the coupling f = when ^ < 0. The simplest choice f = Xcj) makes the perturbations u± 
diverge at t = m/X. However, this coupling is perfectly acceptable in cosmological models 
whose initial singularity is located at f > m/X. For example, if the Universe is born at f = 0, 
the coupling f = X(j) with A < 0 yields an everywhere regular dynamics. 

Before I close this section, let me note that matter-inflaton coupling can signihcantly 
modify the dynamics of matter helds. In particular, the particle production rate can be 
increased. For example, the coupling f = Xcf) with A < 0 diminishes the value of the friction 
coefficient to 


A = 3H + 


X 


m — Xt 

If |A| is large enough, there is a time interval in which the friction becomes negative. During 
that time, the initial matter perturbations keep growing, despite the rapid expansion of the 
Universe. Obviously, this makes the production of particles more effective. Unfortunately, 
the negative friction in this example does not last long. Compared with the typical infla¬ 
tionary period, it is more than one hundred times shorter. Still, it is always possible to 
adjust the matter-inflaton coupling to obtain the desired particle production rate. 


VI. EXAMPLES 

In this section, I shall analyze three specihc choices of the Universe dynamics. In the 
hrst, a toy model is used for the demonstration of how the described procedure works in 
practice. In the second, I propose an inflationary model close to the standard model of the 
Universe. The third is a model of the bouncing Universe. Neither of these models considers 
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matter fields. They are solely used for the demonstration of how a dark energy action is 
associated with a chosen background geometry of the Universe. Hopefully, the right choice 
of the scale factor a{t) will ultimately lead to a realistic cosmological model. 


A. Toy model 


Let me consider a homogeneous, isotropic and spatially flat geometry (fTOj) with the scale 
factor of the simple form 


Its graph is displayed in 
epoch begins at inhnite 


a{t) = 


(40) 


cosh ut 

Fig. O It describes an ever existing Universe, whose inflationary 
past and lasts inhnitely many e-folds. The exit from inflation 



FIG. 2: Toy model. 


happens at t ~ —l/co. The constant a; is a free parameter of the model. 

The scale factor (1401) is a solution of the sigma model (IT^ in which the potential V (0) and 
the target metric F{(j)) are calculated from flTTl) and flT^ . This procedure straightforwardly 
leads to 


F(0) = 


2a;^ 


U(0) = 


6 sinh ud) — 2 


■ u 


(41) 


cosh^ oj(f) ’ cosh^ oj(f) 

With U(0) and U(0) dehned by flTTl) . the action functional flT^ has the solution flTTl) . in 

which gld'J is dehned by ffTOj) and fITOD . The model can further be simplihed by a suitable 
redehnition of the scalar held. Specihcally, the redehnition 0 —)■ x{4>) of fho form 


X = 2 arctan (sinh u(j)) 


(42) 


brings the action flT^ to the form 


R - - Uix) 
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The new potential Uix) = ^(0(x)) reads 


Uix) = 2u^ (4sin2|-l). 

In terms of x, the nontrivial vacuum solution (p = t becomes 

x(t) = 2 arctan (sinha;t), (43) 

while the metric = gjjP) remains the same. The graphs of the potential U{x), and 
the solution (143|l are displayed in Figs. |3] and 01 It is seen that the potential U{x) is a 



FIG. 3: Potential function U{x)- 

periodic function of y, with the period 27r. Thus, the scalar held y lives on a circle. The 
soliton solution (HHll is one-to-one mapping —>■ S^. As for the trivial solutions, the 
theory accommodates two of them. The hrst is given by y = 0, and Anti-de Sitter metric. 
The second has y = ivr, and the metric is de Sitter. Only the hrst solution is stable. 



FIG. 4: Soliton solution. 


because y = 0 is the minimum of the potential Uix), and Anti-de Sitter metric is a stable 
solution of the corresponding geometric equation [38| . The unstable solution y = ± tt is the 
limiting case of the stable soliton solution when t —)■ ± cxd. It seems as if stable soliton is 
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asymptotically unstable. However, there is no contradiction in this unusual situation. This 
is because the stability analysis of Sec. HV] deals with the inhnitesimal perturbations, which 
preserve the monotonous character of the soliton solution. As a consequence, the scalar y 
can be gauged away. This is not the case with the trivial vacuum y = ± vr. No matter 
how small the perturbations of this vacuum are, they can not be gauged away. This is why 
the stability of y = ivr can not be a substitute for the asymptotic stability of the soliton 
vacuum. 

Before I close this subsection, let me say something about perturbations of the background 
solution fHU]) . It is immediately seen that t = 0 is the only critical point of this model. In 
this point, the coefficients of the tensor equation fl23p are regular, but the friction coefficient 
of the scalar equation fl7r|) diverges as —2/t. Luckily, this singularity turns out not to 
be harmful. Indeed, a careful analysis shows that scalar perturbations formed in the past 
regularly pass the critical point t = 0. 

Finally, let me calculate the friction coefficient of the scalar mode of this model. A 
simple calculation shows that it approaches the value 5a; when t —)■ — cxd. In the same epoch, 
the friction coefficient of matter scalars (which obey the equation (133|1 ) has the value 3c;. 
Thus, during rapid expansion of the Universe, the inflaton decays faster than matter scalars. 
Moreover, if the model is dehned by 

_ 

a{t) = (coshcut) " , 

the ratio of the two friction coefficients becomes 1 + 2n/3. When n S> 1, the inflaton decay 
rate becomes much larger than that of matter scalars. 


B. Inflationary Universe 


In the second example, I shall examine the scale factor of the form 


a(t) 


1 + tanh(8a;f) 


In 


+ exp 


cat — 4\ 

4 ;■ 


(44) 


As seen from its graph in Fig. 5, it mimics the standard model of the Universe. Indeed, 
all the expected phases of the cosmological evolution are there. The inflationary epoch 
begins at inhnite past, and lasts inhnitely many e-folds. The exit from inflation happens at 
ut ^ 0, when the early acceleration stops. The Universe continues to expand slowly, until 
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FIG. 5: Inflationary Universe. 


it reaches the moment when the late time acceleration begins. The present epoch is located 
at ojt 7.7. I shall demonstrate later how the parameter oj and the present time tnow are 
calculated from the known values of the Hubble and deceleration parameters. 

The action functional whose vacuum solution is dehned by (1441) has the form (|T^ . with 
F(0) and H(0) calculated from flTTll and flT^ . A straightforward procedure leads to the 
cumbersome expressions which I choose not to display here. Instead, their graphs are pre- 



FIG. 6 : Target metric F((()). 

sented. The action (B, with F[(p) and H(0) depicted in Figs. |6] and [71 has the soliton 
solution (f) = t. 



FIG. 7: Potential V{cj)). 

The parameter u can be determined from the known values of the Hubble and deceleration 
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parameters. These are defined as 


H = 


q = 


a ^ aH"^ 

With the scale factor given by (jS]), the parameters H/uj and q depend on u 
through the combination out. The direct calculation yields the functions whose 


and t only 
graphs are 


-1 


1.5 

1 

0.5 


H/& 

i 2 3 4 5 6 7 S 



FIG. 8: Hubble and deceleration parameters. 

displayed in Fig. | 8 l Using these, one easily finds that the astronomically observed values 

FT = 0.075 Gyr-\ g = -0.5 
imply ut = 7.68 and H/u = O.IA. Thus, 

u = 0.53 Gyr“^, t^ow = 14.55 Gyr. 

The obtained value of fnow is the time coordinate measured from the arbitrarily chosen 
origin f = 0. This is not what one would like to have. Instead, the present epoch should 
be measured relative to a physically significant moment in the history of the Universe. This 
can not be the initial singularity, as our model does not have one. Instead, my choice is the 
end of inflation. The end of inflation t-mi is naturally defined as the moment when the early 
acceleration stops. It is seen from the graph in Fig. [5] that (utinf = 0.003. This leads to 

Gow - Fnf = 14.54 Gyr. (45) 

The time interval flTSl) is a substitute for what is commonly called the age of the Universe. 

So far, I have discussed everywhere regular cosmological models. However, these can 
easily be modified to become models with the initial singularity. For example, the model 
under consideration can be redefined by replacing its scale factor a{t) with 

d{t) = a{t) — a{to). 
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The new scale factor describes a Universe which is born ai t = to, and lives regularly ever 
after. As opposed to the inflationary epoch of everywhere regular model (jS]), the inflationary 
epoch of the new model lasts a hnite number of e-folds. This number can be made arbitrarily 
large by letting to —)■ —oo. It can also be shown that the new model has no singularities 
other than t = to- 

Finally, let me briefly discuss the propagation of small perturbations. It is seen from 
Figs. |5] and [6] that neither the Hubble parameter H, nor the target metric F have zeros. 
As a consequence, the model under consideration has no critical moments. This means that 
both, tensor and scalar, perturbations have everywhere regular dynamics. In particular, 
their friction coefficients are everywhere hnite. While tensor huctuations have always positive 
friction, the friction of the scalar huctuations can become negative. Indeed, a straightforward 
analysis shows that, after the inhation, the scalar friction abruptly drops and becomes 
negative. The period of negative friction does not last long. When the late time acceleration 
begins, it already has a small positive value which gradually approaches zero as t ^ oo. In 
the next subsection, I shall present the example of a model whose scalar huctuations have 
critical moments. 

C. Bouncing Universe 

In this example, I shall consider the scale factor of the form 

a{t) = (l -|- . (46) 

Its graph is displayed in Fig. [9l It dehnes a bouncing Universe which begins in the inhnite 



FIG. 9; Bouncing Universe. 

past, slowly shrinks to its minimal size, and then bounces to an expanding phase. The action 
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functional whose vacuum solution is defined by has the form flT^ . with -F(0) and ld(0) 
calculated from flTT]) and flT^ . A straightforward procedure leads to 


F{cj)) = 


u'^cjF - 1 

( a ;202 + 1)2 > 


V (0) = 


+ 2 

2 (a; 2(/)2 + 1)2 • 


(47) 


The action flT^ . with -F(0) and V{(j)) dehned by flT7)) . has the spatially homogeneous and 
isotropic solution 


(j){t) = t , = —dt'^ + \/l + [dx^ + dy‘^ + dz^') . 


It is seen that the background spacetime is flat in the inhnite past and future. Indeed, all 
the curvature invariants are shown to fall off as 1/F or faster as |t| —)■ oo. Thus, the Universe 
in this example evolves out of the flat spacetime. 

How do metric perturbations propagate in this background? A simple analysis shows that 
tensor perturbations are everywhere regular, whereas scalar ones have critical points. It is 
straightforward to verify that there are three critical points: the two zeroes of id, located 
in t = il/o;, and the zero of id, located in t = 0. These are the singularities of the friction 
coefficient in fl27)) . Its behavior in the vicinity of the critical points is given by 

{ 2 

— in the vicinity of t = 0 , 

/ 

-— in the vicinity of f = =F l/ca. 

t ± l/u 

As has already been mentioned, the singularity of the form —2/t does not violate the regular 
propagation of scalar perturbations. However, it is not the case with the singularities t = 
Fl/uj. These are shown to act as barriers to the scalar perturbations coming from the 
past. Indeed, the scalar perturbations go to inhnity as they approach t = il/cj. This is 
an undesirable property, which calls for the rejection of the model. Still, one should be 
aware of the fact that our linear analysis makes no sense if the helds are too strong. In such 
situations, the interaction terms should be taken into account. Hopefully, this could cure 
the singularity problem. 


VII. RECAPITULATION 

The purpose of this work has been to demonstrate how an arbitrarily chosen background 
of the Universe can be made a solution of a simple model. To this end, I made use of the 
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concept of geometric sigma models. These models possess two distinctive features. The 
first is that any metric can be made a solution of a particular geometric sigma model. The 
second ensures that the complete matter content can be gauged away. In this paper, a 
geometric sigma model is associated with an arbitrary homogeneous, isotropic and spatially 
flat geometry. In its simplest form, the model describes one scalar field in interaction with 
Einstein’s gravity. It possesses the vacuum solution 0 = t, gfj.u = gjiJ- It is important to 
emphasize that, while the background metric gj^J can be chosen arbitrarily, the physics of its 
small perturbations can not. In fact, the role of the background metrics gl^u is to parametrize 
the class of models presented in this paper. This way, the search for a viable cosmological 
model reduces to the proper choice of the background metric. 

The present work begins with the recapitulation of the concept of geometric sigma mod¬ 
els. The construction of the generic model is presented, and subsequently applied to the 
homogeneous, isotropic and spatially flat geometry of the Universe. Then, the dynamics 
of small, localized perturbations of the vacuum is examined. It is demonstrated how all 
but three degrees of freedom can be gauged away. The classical stability of the gauge hxed 
linearized theory is proven in Sec. IIVI This is done by direct calculation, as stability of the 
vacuum solution is not guaranteed by the very construction of geometric sigma models. The 
vacuum stability against matter fluctuations is considered in Sec. |Vl It is shown that the 
inclusion of matter helds does not compromise the results of the preceding sections. The 
analysis of the stability against matter perturbations concludes the general considerations 
of the paper. 

The rest of the paper is devoted to examples. The hrst is a toy model used to demonstrate 
how the method works in practice. The second is the example of an inflationary Universe, 
and the third is a bouncing model. The corresponding action functionals are constructed 
along the lines described in Sec. [TTl The obtained target metric ^(0), and the potential 
V (0) are presented through their graphs. The graphical method is also used for numerical 
calculations. In particular, it is demonstrated how the parameters of the model are calculated 
from the known values of the Hubble and deceleration parameters. 

I have also given a brief insight into the propagation of small perturbations. It is argued 
that the generic linear stability proven in Sec. IIVI fails in the vicinity of critical moments 
H = 0. There, the scalar perturbations diverge, so that a consistent stability analysis must 
go beyond the linear approximation. A more elaborate analysis of this topic is left for the 


future investigation. In particular, the nonperturbative analysis near the critical moments, 
and the study of more complex geometric sigma models is planed. 

At the end, let me point out that the procedure described in this paper misses an im¬ 
portant ingredient. It concerns the back reaction of quantum vacuum fluctuations on the 
background geometry. Indeed, the only matter that geometric sigma models of Sec. [ITl deal 
with is a specihc dark energy of purely geometric origin. Although it successfully generates 
any desirable geometry of the Universe, the influence of quantum fluctuations of ordinary 
matter has not been taken into account. Instead, it has been demonstrated in Sec. El that 
matter helds preserve the classical linear stability established in Sec. IIVI The completion of 
the incomplete cosmology presented in this paper requires the inclusion of quantum fluctu¬ 
ations of both geometry and matter. Until then, it is comforting to know that dark energy 
is commonly believed to dominate all other forms of matter in the Universe. Owing to this, 
the predictions of this work may not be far from realistic, after all. 
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